Solution of the Anderson impurity model via the functional renormalization group 
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We show that the functional renormalization group is a numerically cheap method to obtain 
the low-energy behavior of the Anderson impurity model describing a localized interacting electron 
coupled to a bath of conduction electrons. Our approach uses an external magnetic field as flow 
parameter, partial bosonization of the transverse spin fluctuations, and frequency-independent in- 
teraction vertices which are fixed by Ward identities. We calculate the quasi-particle residue and 
the spin susceptibility in the particle-hole symmetric case and obtain excellent agreement with the 
Bethe ansatz results for arbitrary strengths of the interaction. 

PACS numbers: 72.15.Qm, 71.27.+q, 71.10.Pm 



The Anderson impurity model (AIM) describes a lo- 
calized interacting electron in contact with a bath of 
non-interacting conduction electrons^. The model was 
first introduced in the context of material science for de- 
scribing the emergence of local moments in metals and 
has been studied for half a century by means of a vari- 
ety of methods^. In the past decade, renewed attention 
has been drawn to the AIM because of its experimen- 
tal realization in quantum dot systems in contact with 
a metallic environment. Moreover, the solution of the 
AIM is one of the fundamental steps in the so-called dy- 
namical mean-field theory^. In practice, this step is of- 
ten carried out using Wilson's numerical renormalization 
group (NRG), which yields numerically controlled results 
for the thermodynamics and the spectral propertied. In 
the 1980s the thermodynamics of the AIM has also been 
obtained exactly via the Bethe ansatz^]. For later refer- 
ence, we quote here the Bethe ansatz results for the spin 
and charge susceptibilities in the particle-hole symmetric 
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Here u = U/(irA), where U is the interaction between 
two localized electrons with opposite spin, and A is the 
hybridization energy to the conduction band, in the limit 
where the latter has infinite bandwidth and constant den- 
sity of states. For u > 2 the charge susceptibility \c be- 
comes exponentially small, while the spin susceptibility 
Xs = (4Tk) _1 is proportional to the inverse of the Kondo 
temperature T K = A(Tm/2) 1 / 2 e-^' u / 6+1 ^ 2u \ 

Although Eqs. ( la lb]) can be confirmed numerically 



using the NRG (up to small deviations due to the finite 
width of the conduction band in the NRG calculations), 
it would be useful to have an alternative and numerically 
less expensive method to obtain the correct low-energy 
physics of the AIM at strong coupling, which is domi- 
nated by the exponentially small Kondo temperature. In 
this work we show that this can be achieved by means of 
the functional renormalization group (FRG) methocP. In 



the past few years, many authors have studied the AIM 
using different versions and truncations of the FRGP^, 
but none of these works succeeded in reproducing the 
correct strong coupling behavior of the AIM. We show 
below that this problem can be solved by means of a 
rather simple truncation of the FRG hierarchy involving 
only frequency-independent vertices, provided that the 
transverse spin fluctuations are properly bosonized and 
the interaction vertices are fixed by Ward identities. 

For simplicity, we focus on the particle-hole symmetric 
AIM. After integrating out the conduction electrons, we 
arrive at the following Euclidean action for the Grass- 
mann field d a {r) describing the localized electron, 
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Here T is the temperature (we set kg = h = 1 and focus 
on the zero temperature limit below), = denotes 
summation over fermionic Matsubara frequencies ioj, 
n a = d (J (r)c? cr (r) represents the number of localized elec- 
trons with spin a =t,4o and d ua — J^ T dr e lWT d CT (r). 
Assuming that the conduction band has infinite width 
and constant density of states, the inverse Green function 
is Gq\{vjS) = w — £ CT + iAsgnw, where in the particle- 
hole symmetric case the energy of the localized electron 
relative to the chemical potential is £ CT = —U/2 — aH , 
and H is the external magnetic field in units of energy. 

Our strategy is inspired by the field-theoretical renor- 
malization group approach recently developed by Ed- 
wards and Hewsoi£3 who noticed that for large mag- 
netic fields the AIM can be studied perturbatively, even 
for U ^> A, as long as U/H is sufficiently small. Fol- 
lowing Ref. [T5], let us therefore use the external mag- 
netic field H as a flow parameter for the FRG. Unfortu- 
nately, using this magnetic field cutoff scheme we have 
not been able to reproduce the strong-coupling physics 
of the AIM if we formulate the FRG exclusively in terms 
of fermionic degrees of freedom. The problem is that 
the resulting fermionic vertices exhibit a rather strong 
frequency dependence which cannot be neglected. To 
avoid this problem, we partially bosonize the action |2]) 



by decoupling the interaction in terms of a complex bo- 
son field x representing transverse spin fluctuations. It 
turns out to be advantageous, however, to bosonize only 
a part of the interaction term, while retaining a fermionic 
parametrization for the rest. Therefore, we write the in- 
teraction in Eq. ([2]) as Un^n^ = U\\n\n^,— U_\_s+S-, where 
U — C/|| + U±_ and the Grassmann bi-linears s + — d^d^ 
and S- = d\.d^ represent the usual spin ladder operators. 
Only the transverse part —U±s+S- is now bosonized by 
introducing a complex boson field x via a usual Hubbard- 
Stratonovich transformatiorP^-4 We choose U± such 
that it vanishes for H — > 00, and that it reduces to 
the bare interaction U for H — >• 0. Therefore we set 
[/J 1 = U^ 1 + Rh and require that Rq = and i?oo = 00. 
Note that Rh plays the role of a regulator function for the 
bosonic sector of the theory which switches on the inter- 
action between transverse spin fluctuations. The proper 
choice of Rh will be discussed in more detail below. 

It is now straightforward to write down formally ex- 
act FRG flow equations for t he ir reducible vertices of 
our mixed boson-fermion modePESD For our purpose, it 
is sufficient to neglect the frequency dependence of all 
one-line irreducible vertices with more than two external 
legs. We also neglect the contribution from charge- and 
longitudinal spin fluctuations to the FRG flow^. The 
resulting FRG flow equation for the self-energy Y, a (ibS) 
of the localized electron is shown graphically in Fig. [T] 
the corresponding analytic expression reads 
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The single-scale propagators in our cutoff scheme are 

G a {iu) = -aGl{iui), F ± (iu) = -Fl{iQ)dR H /dH, (4) 

where G a (iu>) is the exact fermionic Green function 
and F±(iu>) is the exact propagator of the spin fluc- 
tuation field x- These propagators can be expressed 
in terms of the corresponding irreducible self-energies 
H a (iuj) and H±{iui) via the Dyson equations G~ 1 (ia;) = 
G^(iu) - and FJ 1 ^) = f/J 1 - U ± (icj). We 

denote bosonic Matsubara frequencies by iu). The right- 
hand side of the flow equation (|3| depends on the ir- 
reducible four-point vertices T 4 " "' "and Y d ° d ' T * x at 
vanishing external frequencies, and on the three-point 
vertices Ig = p^ rf tx _ ^ rp ne SU p erscr ipt s indicate 
the types of external legs attached to the vertices. 

Our goal is to derive flow equations for the spin- 
dependent part M of the fermionic self-energy at van- 
ishing frequency and the quasi-particle residue Z, which 
are defined via the low-energy expansion 
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FIG. 1. (Color online) Graphical representation of the FRG 
flow equation |3| for S^(ioj) (shaded circle). Solid and wavy 
lines represent full fermion and boson propagators, respec- 
tively. For each type of line, the slash insertions denote the 
corresponding single-scale propagators in Eq. Q. The dot on 
the left-hand side represents a derivative with respect to H. 



In order to obtain M and Z from the solution of the 
flow equation ^ we need to know the bosonic self-energy 

U±(iQ) and three different vertices: rf^""**"**, rf d " xx , 
and 7. One could write down additional flow equations 
for these vertices, but these depend on higher order ver- 
tices and it is not clear how to truncate this infinite hier- 
archy. We shall therefore adopt a different strategy which 
uses Ward identities and skeleton equations to obtain a 
closed system of flow equations for M and Z. 

First of all, we note that the three-point vertex at van- 
ishing frequencies can be related to the self-energy pa- 
rameter M via the Ward identity 



7 = b/(h + u± arctan b) , (6) 



which can be derived for our bosoni zed t heory follow- 
ing the work by Koyama and Tachikf^EU. Here uj_ = 
U±/(nA) and b = h + m, with m = M/A and h = H/A. 

Next, let us derive an alternative flow equation for 
M which does not explicitly involve the four-point ver- 
tices. To this end we consider the longitudinal spin- 
susceptibility x\\ — ds/dH, where s = (n^ — n±) /2 is 
the local magnetic moment (which we normalize such 
that s = 1/2 corresponds to a fully magnetized state). 
The local moment is then connected to the self-energy 
parameter M via the Friedel sum rule^, 



its = arctan(/i + m) = arctan b. 



(7) 



Taking the derivative with respect to h we obtain 

dm/dh = x\\/p-l, (8) 

where p = [7rA(l + & 2 )] -1 is the exact density of states 
at vanishing energy. To determine x\\ we use t ne Ward- 
identity^! 

X\\/p = l/Z + P r ± , (9) 

where T± is the exact interaction vertex between two 
electrons with opposite spin at vanishing frequencies. In 
our partially bosonized theory we have 



£ CT M = U/2 - oM + (1 - Z- 1 )^ + 0{lo 2 ). (5) 
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But the three-point vertex and the fermionic four-point 
vertex turn out to be mutually related via a skeleton 
equation. With the help of standard functional meth- 
ods^ it is indeed straightforward to show that, at finite 
frequencies, 

Tp dlX (uj + ui,uj;ui) = 1 - [ G^')G t (iu/ + iu>) 

J UJ' 

xrf dM {cj + u, u, oj' + oj). (11) 

Using the above relation to express p^ d 4A d t m terms f 
Y<^t d ix = ^ a ^ vanishing external frequencies, we find for 
the longitudinal susceptibility^, 



P 



1 

Z 



1 



1 + b 2 arctan b 



(12) 



where we have used the Ward identity ( |6[) to eliminate the 
vertex 7. Substituting the expression ([12]) into Eq. pL 
we obtain a flow equation for m which depends only on 
the self-energy parameters m and Z. 

Finally, to close our system of flow equations we ap- 
proximate the bosonic self-energy Tl±(iuj), appearing im- 
plicitly in Eq. (jsj), by 

nAn ± {iu>) = *fP ± (\Q\/(ZA),bagau), 



where the dimensionless function 

ln[l + x/(l — ib)] 2i arctan b 
{x -2ib)(l + x/2-ib) x - 2ib 



P±(x,b) 



(13) 



(14) 



is obtained from the skeleton equation^ for H±(iu>). The 
prefactor in Eq. (13) is determ ined by demanding that 
n_i_(0) satisfy the Ward identity^!!!! 



n ± (fi)/[i-u±u ± (o)] = x± = 8/H 



(15) 



relating the transverse spin-susceptibility x± to the local 
moment s. Using the Friedel sum rule Q it is easy to 
see that our expression (l2"| for the longitudinal suscep- 
tibility is consistent with Eq. (15) for \Xi m the sense 
that for H — > both susceptibilities approach the same 
limit Xs- We point out that Eqs. §6§ and ( 15 1 imply that 
a = [1 — U±U±(0)} = h/(h + u_LarctanF) > 0, which 
guarantees that there is no magnetic instability. 

We have now obtained a closed system of flow equa- 
tions for m = M /A and Z. For simplicity, we expand the 
function P±(x,b) in Eq. (13) to linear order in x, which 
is consistent with the low-energy expansion ([5| in the 
fermionic sector. After straightforward algebra we arrive 
at the flow equation hdZj dh = rjZ for the quasi-particle 
residue Z ' . with 
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Here r = ttARh is the dimensionless regulator function 
of the boson propagator, and 

c = (u_ L 7)(2^)" 1 [(1 - ib)~ 2 - arctan b/b] . (17) 




FIG. 2. (Color online) Logarithmic plot of the bosonic regu- 
lator r = ttARh as a function of the magnetic field h = H/A, 
for U/{irA) = 1 (solid line) and U/(nA) = 8 (dashed line). 
The regulator is linear in h both at small and large h, with 
a change in the prefactor from a (h < 1) to /3 (h > 1). The 
inset shows the [/-dependence of the coefficients a and j3. 



The above system of flow equations can be easily solved 
numerically once the bosonic regulator function r = 
ttARh has been specified. The simplest choice of a lin- 
ear magnetic field dependence of the bosonic regulator 
(r cx h, in agreement with the magnetic field regulator 
in the fermionic sector) does not lead to satisfactory re- 
sults. After testing many regulators, we found that the 
best results can be obtained with a regulator of the form 



r= (a-/3)h/(l + h 2 )+ph, 



(18) 



which is linear in h both for small and large h, but ex- 
hibits a change in slope from a to at h « 1, as shown 
in Fig. [2j Comparing our FRG results at weak coupling, 
U A, with the known perturbative expansion for Z 
and Xs , we also realized that the coefficients a and /3 are 
weakly dependent functions of the dimensionless bare in- 
teraction u — U/(ttA). Their functional dependence is 
well described by rational functions of low degree. Hence, 
we introduce the following ansatz, 

a = a /{a 2 + u 2 ), /3 = p x u/{p 2 + u 2 ), (19) 

where the four parameters ao, 0.2, fix, and P2 can be 
determined by matching the FRG results for Z and Xs at 
weak coupling (e.g., for u = 0.1) with the corresponding 
values obtained in the perturbative expansion at 0(u 2 ), 
and by imposing that the Wilson ratio R is equal to 2 
for two values of u in the Kondo regime (e.g., we impose 
R = 2 for u = 4 and u = 8). Note that in the AIM the 
expression for the Wilson ratio reads 



p _ X s(U)/Xs(0) a 7 

R =z=mi^m = ^ XsZ= x. 



2x s 
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where the last equality follows from the Yamada-Yosida 
Ward identitie a 19 * 21 !. At large U, charge fluctuations are 
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FIG. 3. (Color online) Inverse quasi-particle residue Z~ x and 
spin susceptibility \a as functions of U / (ttA) for H — 0. The 
FRG results for Z -1 (squares) and \s (circles) are compared 
with the corresponding Bethe ansatz results (solid and dashed 
lines, respectively). As a reference, we also show Z~ x calcu- 
lated within second order perturbation theory (dotted line). 
The inset shows the FRG (diamonds) and Bethe ansatz (solid 
line) results for the Wilson ratio R. 



strongly suppressed and the low-energy physics of the 
AIM is effectively the same as in the Kondo model. In 
this regime the charge susceptibility Xc is therefore negli- 
gible and the Wilson ratio takes the Kondo-model value 
R = 2. For the parametrization of a and /3 in Eq. (191, 
we find cxq 



0.36, a 2 « 1.4, ft « 0.19, and (3 2 « 5.4. 
The Independence of a and (3 is shown in the inset of 
Fig. U 

Having fixed the regulator in the bosonic sector, we 
can now solve our FRG flow equations numerically, with 
a minimal computational cost. Choosing the initial mag- 
netic field Hq sufficiently large, we may use as initial con- 
ditions for M and Z the Hartree-Fock values Mo = U/2 
and Zq = 1, which are exact in the limit Hq —¥ oo. The 
FRG results for Z and \s > in the limit of vanishing mag- 
netic field, are shown in Fig.[3j where in the inset we also 
show the corresponding Wilson ratio. The remarkable 
agreement of our approach with the Bethe ansatz results 
for arbitrary large values of the interaction shows that the 
present FRG truncation, supplied with exact Ward iden- 
tities, is able to capture the exponentially small Kondo 
scale in the low-energy quasi-particle resonance. 

In Fig. [4] we show, for u — U/(ttA) — 5, the magnetic 
field dependence of Z, s, and 7. In the magnetization 
curve s(H) one can clearly identify the two energy scales 
characterizing the strong coupling regime of the AIM, 
namely the bare Coulomb interaction U and the width 
ZA of the Kondo resonance. Indeed, reducing the exter- 
nal magnetic field H from a large initial value, at H ps U 
we start observing a slow decrease in the magnetic mo- 
ment from the fully magnetized value s — 1/2. The de- 
crease in s becomes then faster when H sa ZA, until the 
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FIG. 4. (Color online) Inverse quasi-particle residue Z~ x (dot- 
dashed line, scale on the right), local magnetic moment s 
(solid line), and three-point vertex 7 (dashed line) as functions 
of the external magnetic field H, for U/(nA) — 5. 



magnetization eventually vanishes for H -C ZA. The 
presence of two different energy scales is also reflected 
in the non-monotonic behavior of the three-point vertex 
7(i?). Starting from the bare value 7 = 1 at large H, the 
vertex becomes initially larger than unity for H U and 
later decreases again, approaching (in the Kondo regime) 
the value l/u < 1 for H < ZA. 

In summary, we have developed an FRG approach to 
the AIM which uses the external magnetic field H as a 
physical flow parameter, and where the fermionic inter- 
action U is partially decoupled via a bosonic field de- 
scribing the transverse spin fluctuations. We have trun- 
cated the exact hierarchy of FRG flow equations keeping 
only the frequency- independent parts of the interaction 
vertices and expanding the fermionic and bosonic self- 
energies to linear order in frequency. With an extensive 
use of Ward identities, we have been able to express all 
the relevant vertices in terms of self-energy parameters, 
thereby avoiding further approximations in the flow equa- 
tions for the vertices. Comparing our results with the 
exact Bethe ansatz solution at vanishing magnetic field, 
we have shown that the present FRG approach is able to 
reproduce the exponential [/-dependence of the width of 
the Kondo resonance. Moreover, the use of the magnetic 
field if as a physical cutoff for the FRG flow equation 
gives access to the non-trivial ii-dependence of physical 
observables such as the local magnetic moment. Possi- 
ble extensions of the present method include the study 
of the non-symmetric AIM, the use of finite frequency 
terms in the fermionic self-energy in order to describe 
the high-energy features of the spectral function, and a 
generalization of our approach to non-equilibrium situ- 
ations relevant for the transport properties of quantum 
dots. 
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